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Abstract 

£^ ■ A sensitivity of twist-4 and a s values extracted in the NLO QCD analysis of nons- 

inglet SLAC-BCDMS-NMC deep inelastic scattering data to the choice of QCD renor- 
malization scale (RS) is analysed. It is obtained that the high twist (HT) contribution 
to structure function F2, is retuned with the change of RS. This retuning depends on 
the choice of starting QCD evolution point Qq and x. At Qq > 10 GeV 2 the HT con- 
Q I tribution to F2 is retuned at small x and almost not retuned at large x; at small Qq it 

■ exhibits approximate RS stability for all x in question. The HT contribution to is 

RS stable for all Qq and x. The RS sensitivity of a s also depends on the choice of Qq: 
at large Qq this sensitivity is weaker, than at small one. For Qq = 9 GeV 2 the value 
a s (M z ) = 0.1183 ± 0.0021(stat + syst) ± 0.0013(RS) is obtained. Connection with the 
higher order QCD corrections is discussed. 
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1. Interest to the quantitative description of high twist (HT) contribution to the deep 
inelastic scattering (DIS) cross sections has been increased recently, in particular, due to 
the development of infrared renormalon (IRR) model (see e.g. review |]J). Within this 
model one can derive the x-shape of HT contribution from the x-shape of leading twist 
(LT) structure functions. This connection allows for to obtain precise predictions for the HT 
contribution since the LT contribution can be rather precisely determined from experimental 
data. Meanwhile the experimental determination of HT contribution is not direct and is 
based on fitting a combination of log- and power-like terms to the data. If the data accuracy 
is not high enough, the correlation between these log- and power-like terms can be large, that 
was explicitly shown in the combined SLAC-BCDMS data analysis of Ref. 0. If the large 
correlations occur, the separation of terms is unstable with respect to the various inputs of 
fit and thus it becomes important to study the stability of HT separation. One of the poorly 
defined ansatz of a QCD analysis of DIS data is the choice of renormalization scale (RS). The 
uncertainty due to RS variation is connected with the account of higher order (HO) QCD 
corrections since in the analysis with complete account of HO terms, the RS dependence 
of fitted parameters should vanish and thus the observed RS dependence can be used for 
the estimation of HO terms effect. Earlier we reported the results of NLO QCD analysis of 
high x SLAC-BCDMS-NMC data || . A short communication on the RS dependence of HT 
contribution and a s obtained in this analysis was reported in Ref. |3J]. In this paper more 
detailed study of this dependence is given. 

2. Our approach used for the study of RS stability of DGLAP evolution equation in 
NLO QCD is the same as described in Refs. || |J. Within this approach the RS of QCD 
evolution is changed from Q to k^Q, where &r is arbitrary parameter, conventionally varied 
from 1/2 to 2. This approach contains certain simplification since the change of scale can 
depend on x, but in our analysis this effect is not so essential due to the limited range of 
x. For the nonsinglet case NLO DGLAP equation with an arbitrary choice of RS looks as 
follows: 

Q dqNS ^ Q) = ^iMl P NS,(0) g, g NS + ^RQ) [p NS,(l) g NS + ln(fcR)A) pNS,(0) ^ g N S] 

(i) 

where P®q = f dzP(z)q(x/z, Q) denotes Mellin convolution; the evolved distribution; 

pNS,(o) an( j pNS,(i) are lo and NLO parts of splitting function; a s (Q) is the running 
strong coupling constant; and /3 is the regular coefficient of renormgroup equation for a s : 

n da s _ _Po_ 2 _ 3 

The equation ([!]) was solved with the help of direct integration method implemented in the 
code used earlier ||. In the analysis of Ref. |J we used combined SLAC-BCDMS-NMC 
proton-deuterium data @ with the cuts x > 0.3 to reduce QCD evolution to the nonsinglet 
case and x < 0.75 to reject the region where a nuclear effects in deuterium can be significant. 
The initial scale of evolution was chosen equal to = 9 GeV 2 to provide comparability 
with the earlier results of Ref. ||. A complete account of point-to-point correlations due 
to systematic errors was made through the covariance matrix approach, similarly to our 
earlier analysis of Ref. f|. The formula were fitted to the cross section data to allow for 
simultaneous and unbiased determination of a twist-4 contribution to structure functions 
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Figure 1: The values of proton and deuterium H 2 (x) for different values of &r (open circles: 
/cr = 1/2; full circles: k R = 1. 
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Figure 2: The dependence of proton Ag NS on a 2 at = 9 GeV 2 (full lines: log 2 k-& 
dashed lines: log 2 k R = —1/2). 
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and F 2 : 

(i - v - ^) F 2 ^ Q)+ (l- 2 |) f./--,, Q) , 

2xFf r {x, Q) = if T (x, Q) - F L HT (x, Q), 

F 2 H L T (x, Q) = F™ c (x, Q) + ^ 2iL (x)l^, 

where F™ c are the LT contributions obtained as a result of integration of Eqn. ([!]) with the 
account of target mass corrections ||; s is the total c.m.s. energy; m\ is the scattered lepton 
mass; and y is the regular lepton scattering variable. The values of functions H 2 ^(x) at x = 
0.3, 0.4, 0, 5, 0.6, 0.7, 0.8 were fitted, between these points H 2 ,\ J {x) were linearly interpolated. 
The functions H 2 for proton and deuterium were fitted independently, while the functions 
Hl for proton and deuterium due to limited accuracy of data, turned out to be compatible 
within errors and all fits were performed under constraint H^(x) = H^(x). 

3. The proton and deuterium H 2 (x) for &r = 1 and k-& = 1/2 are given in Fig. [I]. One 
can see that they depend on fc R at x ~ 0.3 and practically does not depend at x ~ 0.8. 
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Figure 3: The dependence of Q 2 [a 2 (Q) — a 2 (Qo)] on Q 2 (full line: a s (Qq) = 0; dashed 
line: a s (Qo) = 0.1; dotted line: a s (Qq) = 0.2). The dashed-dotted line corresponds to 
Q 2 [a s 3 (Q) - a s 3 (Q )] for a s (Q ) = 0.18. 

To give explanation of this behaviour recall the basic properties of solutions to the DGLAP 
evolution equations. After linearization on \nk R Eqn. ([!]) can be analytically solved in the 
Mellin momentum space: 

M NS (n, Q) = M NS (n, Q )M NLO (n, a) exp [g(n) (a 2 - a 2 ) ln(Afe)] , 

where a = a s (Q), ao = a s (Qq); M NS (n,Q) are Mellin moments of g NS ; M^o{n,a) defines 
NLO evolution of these moments; g(n) is a linear function of Mellin moments of the splitting 
functions Pqq'^ and P<^ s '^. Introduce a function Aq NS (k R ) = g NS (/cR,) — ? ns (A;r = 1) that 
is convenient to study RS dependence. The Mellin moments of this function M A (n,Q) are 
given by 

M A (n,Q) = M NS (n,Q )M NLO (n,«) {exp [g(n) (a 2 - ajj) ln(k R )] - 1} . (2) 

In Fig. the precise dependence of proton Ag NS on a 2 (Q) for different k R , obtained as the 
result of numerical integration of Eqn. ([[]) at Qq = 9 GeV 2 , is givenF]. The range of a s in 
the figure correspond to the variation of Q 2 from 1 to 9 GeV 2 , i.e. the region where the 
HT contribution is most significant. It is evident that for all values of x in question Ag NS is 
approximately proportional to In k R . This means that the linearization of Eqn. (|1|) is justified 
and that in Eqn. (fj) the part of exponent containing In k R can be expanded, so that 

Ag NS (a;, Q) w In(Aft) {a 2 - a 2 ) q NLO (x, a), (3) 

where qNLo{ x i a ) is the Mellin inverse of the product M NS (n, Qo)M NLO (n, a)g(n). 

The Q-behaviour of ^nlo is defined by M NLO (n, a). At x — 0.3 the function ^nlo weakly 
depends on a (see Fig. 0). This is consequence of the well known effect, that the non-singlet 
QCD evolution has stationary point at x ~ 0.1 due to the fermion conservation. In the 
vicinity of stationary point scaling violation is small, but, due to at large n the function 
A^NLo( n > a) rises with a faster, than at low ones, at large x the scaling violation is more 
pronounced and the function Ag NS rises with a significantly faster than a 2 . 

^^Here and below we do not give the results for deuterium since they are similar to the proton ones. 
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Figure 4: The value of a s at starting evolution point that provides the power-behaviour 
simulation of factors \a^{Q 2 ) — a™(Qo)] in the region from Qi to Q 2 (full line: n — 1, dashed 
line: n = 2, dotted line: n — 3). The value of a s (Q 2 ) is also given for comparison (dashed- 
dotted line). 
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Q^=50 GeV 2 



10 



-2 



10 



10 



-4 




10 

Q 2 [GeV 2 ] 



Figure 5: The dependence of | Q 2 AF| I | [GeV 2 ] on Q 2 at &r = 1/2 and at different values of 
x and Qq (full lines: a; = 0.3; dashed lines: x = 0.8). 



In the NLO non-singlet approximation the LT contribution to F 2 is given by 

2tt' 



F 2 LT = g NS + % NS ®Cf' (1) , 



where C^ IS ' ( ' 1 ^ is the NLO coefficient function. Since the second term of above expression 
is suppressed with respect to the first one at moderate x, the Q-behaviour of function 
AF2(kn) = F 2 lt (A;r) — F^^/cr = 1) at x — 0.3 approximately coincides with the ones 
of Ag NS . The Q-behaviour of the factor [a 2 (Q 2 ) — a 2 (Qo)], coming to Eqn. (H) depends on 
Qq. If Q Qo this factor is ~ l/ln 2 Q, i.e. falls with Q slower, than 1/Q 2 . Meanwhile in 
the vicinity of Qq, where this factor vanishes, its Q-dependence is steeper and it can simulate 
1/Q 2 behaviour in a rather wide range of Q (see Fig. |3]). One can easily show that the value 
of a s at a starting evolution point Qq , that provides the 1/Q 2 behaviour simulation of the 
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Figure 6: The values of H\{x) for different choices of Q and RS (open circles: k R 
full circles: k R = 1; squares: &r = 2. 
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Figure 7: The dependence of Aaf (A;r) = a s (Mz 



-a s (Mz) |fc H =i on the choice of 



renormalization scale for different Qo (open circles: Q$ = 50 GeV 2 , full circles: Ql = 
1 GeV 2 ). For comparison are also given the Ao^ s (/cr) values obtained in the fits at Qq = 
50 GeV 2 with HT fixed at the values obtained in fits with &r = 1 (squares). Error bars are 
not given. 



factor [a™(Q 2 ) — a"(Q )] m the region from Qi to Q 2 , is 
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The dependence of a s [Q^ j on Q 2 for Qi — 1 GeV z obtained with the help of Eqn. (f4|) 
is given in Fig. f|. For the Q 2 interval from 1 to ~ 10 GeV 2 , where the HT are mostly 



Q ( n) 



40 GeV 2 . Due to the weak 



significant, a s yQ j ~ 0.2, that corresponds to 
dependence of ^nlo on Q at small x, the function AF2 also simulates 1/Q 2 behaviour at 
x = 0.3, if the starting evolution point is chosen equal to Q . At fixed Qo this simulation, 
in general, becomes less probable with the rise of x just due to the steeper rise of Ag NS with 
a at large x. 

The behaviour of Q 2 AF% at &r = 1/2 for different x and Qq is given in Fig. [|. At 
50 GeV 2 and x = 0.3 the function AF 2 simulates the 1/Q 2 behaviour about perfectly 
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in the total Q region relevant for HT determination. This leads to the H 2 retuning in the fits 
with different &r, since AF 2 is compensated by the additional contribution to H 2 with the 
sign opposite to AF 2 (see Fig. g). At x — 0.8, due to the fall of <?nlo with Q, the simulation 
is much worse. If Q 2 , = 1 GeV 2 , the absolute value of Q 2 AF 2 steeply rises at x = 0.3 and 
small Q due to the factor [a 2 (Q 2 ) — a 2 (Q )}. At x = 0.8 this rise is not so steep because 
of the fall of ^nlo with Q, but it cannot suppress the general rise. As a consequence, at 
small Qo the function AF 2 cannot simulate 1/Q 2 behaviour at all x and it should be at least 
partially compensated by the change of Q-dependence of the LT contribution. The remnant 
HT retuning is still possible if the Q-dependence of AF 2 at some x is more similar to 1/Q 2 , 
than to the Q-dependence of dF 2 /da s (Mz)f\. The explicit tracing of the balance between the 
H 2 retuning and the change of LT contribution is not so simple due to a s (Mz) is determined 
by data at all x and Q. Anyway, as one can see from Fig. |B| at small Qo the function H 2 
is retuned with the change of /cr significantly smaller than at large Qo- At Qo = 9 GeV 2 
the retuning effect is almost the same as for Qq = 50 GeV 2 . This is natural since the AF 2 
behaviour at small Q weakly depends on Q , if Q is large, and the data are not very sensitive 
to the HT contribution at Q 2 > 10 GeV 2 . Thus the choice of Qo have small effect on the 
fitted HT contribution if Q 2 > 10 GeV 2 . 

The RS stability of a s also depends on the choice of Qo- I n the fit with HT fixed AF 2 is 
compensated by the change of LT contribution, that leads to the shift of a s . If H 2 is released 
in the fit, the shift of a s can change due to the partial compensation of AF 2 by the change 
of H 2 . Since the 1/Q 2 simulation of AF 2 depends on Q , the a s dependence on /cr changes 
with Qo as well as H 2 (see Fig. 0). At large Qo the AF 2 slope on /cr is negative at small 
Q. As a result, in the fit with HT released the H 2 slope on /cr is positive. Correspondingly 
the Q-dependence of LT contribution after HT releasing becomes weaker at large /cr and 
steeper at small ones, i.e. the slope of fitted a s (M z ) value on fc R decreases as compared to 
the fit with HT fixed. In the fit with HT fixed the a s (Mz) slope on /cr is positive and thus 
the RS uncertainty on a s at large Qo becomes smaller. At small Qo the H 2 slope on /cr at 
x ~ 0.5 -T- 0.7 is negative. The data from this region of x have the largest impact on the a s 
determination and thus the H 2 releasing leads to the increase of RS error on a s , although 
the scale of effect is smaller as compared with the fit at large Qo due to the weakness of HT 
retuning at small Qo- In our analysis 

a s (M z ) = 0.1151 ± 0.0015(stat + syst) ± 0.0045(RS) 

for Q 2 = 1 GeV 2 and 

a s (M z ) = 0.1183 ± 0.0021(stat + syst) ± 0.0013(RS) 

for Qq = 9 GeV 2 , where the RS error is estimated as the half of a s (M z ) spread with the 
change of /cr from 1/2 to 2; the central value is shifted to the centre of this spread. The 
values of \ 2 are approximately the same for different /cr, but in view of that at small Qo 
the total a s (Mz) error is about two times larger, than at large one, we consider the a s value 
determined from the fit with large Qo as more reliable. 

The function Ag NS , by definition, is connected with the NNLO QCD corrections to 
evolved distributions. A natural assumption is that the exponent in NNLO part of mo- 
ment expression can be expanded, similarly to the exponent in Eqn. (|2|), and the NNLO 

2 The Q-dependence of LT contribution is mainly driven by a s and thus the balance between the absorbtion 
of AF2 into the LT contribution and H2 is defined by dF2/da s (Mz). 
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Figure 8: The comparison of Q 2 AF 2 p (k R = 1/2) and Q 2 AF[(k R = 1/2) dependence on Q 2 . 

contribution is ~ [a 2 (Q 2 ) — a 2 (Qo)] as well as Ag NS . One can see from Fig. ||| that the 
factor [a^(Q 2 ) — a^(Q )] coming to N 3 LO contribution also can simulate l/Q 2 behaviour. 
Moreover, the region of Q where the simulation is possible widens from the lower orders to 
the higher ones. This gives an indication that the retuning of HT contribution after the 
accounting of HO QCD corrections to DGLAP kernel would exhibit the same properties as 
the RS retuning in NLO. Note that if Q 2 Qi, then a s (Q n ^ ~ a s (Q2) for all n (see 

Eqn. H and Fig. ^). This means, that in the analysis with simultaneous determination of 
a s and HT at large Qo the contribution to F 2 due to the HO QCD corrections to DGLAP 
kernel can be merely fitted together with the HT, if g(n) and the Mellin moments of HO 
splitting functions have the similar large n asymptotes. The fitting of HO corrections sure 
leads to the increase of a s errorQ. Meanwhile the reduction of RS uncertainty, that is one of 
the dominant sources of a s error, is larger and, as one can see from the above, the total a s 
error becomes smaller. 

One can see from Fig. [7|, that at Ql = 50 GeV 2 the fitted a s (M z ) value is a nonlinear 
function of lnA;R contrary to the fits at Ql = 1 GeV 2 . The reason of this difference is the 
large correlation between a s and H 2 (c.f. Ref. @), that depends on how well dF 2 /da s (Mz) 
can simulate l/Q 2 behaviour. With the rise of Qo this correlation increases. For example the 
correlation coefficient p .5 fc> r ^ s (Mz) and H 2 (x = 0.5) at /cr = 1/2 is -0.82 for Ql = 1 GeV 2 
and —0.97 for Ql = 50 GeV 2 . As a consequence, at large Q a small nonlinearity of g NS on 
In k-R better manifests and takes non-negligible effect on the fitted parameters values (remind 
that the effective amplification of nonlinear effects in a fit is proportional to 1/(1 — p 2 )). For 
the comparison, the fits with HT fixed exhibit almost linear dependence of a s (Mz) on ln^R 
(see Fig. |7|). One of the reflections of this nonlinearity is that for Q^ = 50 GeV 2 the difference 
between H 2 (x) at k R = 2 and at /cr = 1 is small (see Fig. ^|). This is unpleasant feature of the 
analysis since the variation range of /cr is conventional and the nonlinearity does not allow 
to rescale the RS uncertainty. One of the possible ways to suppress the nonlinear effects is 
to decrease the correlation between HT contribution and a s , e.g. adding more data to the 
analysis. From another side this correlation leads to the reducing of RS error on a s and it is 
necessary to keep a balance between the linearity and the size of RS error. 

3 In the fit at Ql = 50 GeV 2 with HT fixed the a s (M z ) error is 0.0004 as compared with 0.0021 in the fit 
with HT released. 
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Figure 9: The values of H^(x) for different choices of Qo and RS (open circles: k R = 1/2; 
full circles: k R = 1; squares: k R = 2). 



In NLO QCD the LT contribution to structure function is proportional to the Mellin 

NS,(1) 
L : 



convolution of g NS with the NLO coefficient function C- 



a s NS 

2n 



NS,(1) 



Due to the convolution smearing a function AF L (k R ) = F^ T (k R ) - F^ T (k R = 1) depends on Q 
steeper, than AF 2 . One can see, that at x = 0.3 and large Qo the function Q 2 AF L (k R = 1/2) 
falls in about two times in the region of Q 2 = 1 4- 3 GeV 2 (see Fig. ^). At large x and low 
Qo it falls even more steeply. As a result, H L (x) exhibits only weak dependence on k R for 
all x and Qo (see Fig. ||). 

4. In summary, we can conclude that the HT contribution to structure function F 2 
extracted in the NLO QCD analysis of nonsinglet SLAC-BCDMS-NMC data is retuned with 
the RS change. This retuning depends on the choice of starting evolution point Qo and x. 
At Qo > 10 GeV 2 the HT contribution to F 2 is retuned at small x and not almost retuned 
at large x; at small Qo it exhibits approximate RS stability for all x in question. The RS 
sensitivity of a s also depends on the choice of Qo- at large Qo this sensitivity is weaker, than 
at small ones. The HT contribution to F^ is RS stable for all Qo and x. 

The RS stability of HT contribution is important for clarification of their nature: due 
to the both HT and HO corrections fall with Q, it was often claimed that the extracted 
HT terms can contain contribution from HO. The HT absorbtion by NNLO correction was 
observed in the analysis of neutrino structure function xF% [|l(|, although the effect was 
smashed due to low accuracy of the data. Our results indicate, that in the analysis of high 
statistical charged leptons DIS data the unambiguous separation of twist-4 contribution 
and NNLO QCD corrections to DGLAP kernel is possible if Qo is low. This conclusion 
is especially important because no complete NNLO calculation of the splitting functions is 
available up to now and it is impossible to perform exact direct clarification of this point. 
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